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Stage-2 MnCl2 GIC magnetically behaves like a quasi 2D XY antiferromagnet on the triangular 
lattice (in-plane lattice constant of Mn layer is 3.692 A and the c axis repeat distance is 12.80 A). The 
temperature (T) dependence of the zero-field cooled (ZFC) and field-cooled (FC) magnetizations 
shows a typical spin glass behavior below the spin freezing temperature Tsa (= 1-1 K)- The AC 
magnetic susceptibility shows a peak at T = Tsa at H — 0. This peak shifts to the low temperature 
side with increasing H, according to the de Almeida- Thouless critical line. We have undertaken 
an extensive study on the T dependence of specific heat in the absence of an external field H and 
in the presence of H along the c plane. The magnetic specific heat Cmag at iif = shows no 
anomaly at Tsa, but exhibits double broad peaks around 5 K and 41 K. The anomaly at 41 K is 
the onset of short range spin correlation. The anisotropy of the DC magnetic susceptibility starts 
to become appreciable below 50 K. The magnetic specific heat Cmag at = is described by 
Cmag oc (l/T^) exp(-A/r) with A = 1.41 ± 0.03 K, while Cmag at = 10 kOe is proportional to 
T. The entropy due to the broad peak around 5 K is 1/3 of the total entropy. The residual entropy 
is 63 % of the total entropy because of highly frustrated nature of the system. The magnetic neutron 
scattering indicates that a short range spin order associated with the incommensurate wave vector 
IQincomml i~ 0.522A~^ at 0.45 K) appears below 5 K and remains unchanged down to 63 mK. The 
in-plane spin correlation length is only 18 A at 0.45 K. The low temperature phase below Tsa is a 
kind of reentrant spin glass phase where the spin glass phase coexists with a short range spin order 
associated with \Qincomm\- 



PACS numbers: 75.40.Cx, 75.50.Lk, 75.50.Ee 

I. INTRODUCTION 

In stage-2 MnCl2 GIC, the MnCl2 layer consists of a 
three-layer sandwich of Cl-Mn-Cl layers, with the same 
layered structure as pristine MnCl2. The MnCl2 lay- 
ers are separated periodically by two graphite layers in 
stacks along the c axis. The c-axis repeat distance is 
12.804. The MnCl2 layer forms a triangular lattice (lat- 
tice constant 3.692 A) translationally incommensurate 
with the graphite host but rotationally locked into 30° 
with respect to it. Stage-2 MnCl2 GIC is a suitable proto- 
type for studying the classical two-dimensional (2D) XY- 
like antiferromagnet on a triangular lattice (AFT). The 
interplanar interaction between adjacent MnCl2 layers, 
J' can be greatly reduced by these intervening graphite 
layers, while the antiferromagnetic intraplanar exchange 

interaction between Mn^"*" spins, Ji, remains virtually 
unchanged.EM 3 | 4 | 5 | B | 7 | « i miu | 

The 2D XY-like AFT has received attention from the- 
orists because the spins in it are fully frustrated. The 
ground state consists of spins on three sublattices forming 
120° angles with respect to each other, the x \/3 spin 
structure. Because there are two senses to the spin helic- 
ity, the ground state has a twofold discrete degeneracy as 
well as an XY-like continuous degeneracy. Consequently, 
it is predicted to undergo two phase transitions, one asso- 
ciated with Ising-type symmetry breaking and the other 



with a Kosterhtz-Thouless (KT) mechanism! " ' ^^ ' ^^ 'Note 
that our system has an incommensurate spin structure 
close to 2-\/3 X 2\/3 spin structure (see later). In this 
sence our system is different from an ideal 2D XY anti- 
ferromagnet on the triangular lattice. 

The magnetic properties of stage-2 MnCl2 GIC 
have been extensively studied by DC susceptibility,!^ 
the zero-field cooled (ZFC) and field-cooled (FC) 
magnetization,^ AC magnetic susceptibility,!^ electron 
spin resonance (ESRjP^ specific heat,!^ and magnetic 
neutron scattering .fsEGS! The low temperature phase be- 
low TsQ (= 1.1 K) is a spin glass phase, reflecting the 
fully frustrated nature of the system. A short range spin 
order related to the incommensurate spin structure (typi- 
cally ISA in-plane spin correlation length) appears below 
5 K, and remains unchanged even below Tsg- 

In the present paper, we report the result of specific 
heat measurement of stage-2 MnCl2 GIC between 0.45 K 
and 50 K, with and without an external magnetic field ap- 
plied along the c plane. We show that the magnetic spe- 
cific heat shows no anomaly at Tsg , but exhibits a broad 
peak around 5 K. It also shows a broad peak around 
41 K, which is the onset of short range spin correlation. 
There is a relatively large residual entropy at T = K 
reflecting the fully frustrated nature of the system. The 
results of spcciflc heat measurements will be discussed in 
light of the experimental results of the ZFC and FC mag- 
netization, AC magnetic susceptibility, magnetic neutron 
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scattering, ESR, and DC magnetic susceptibility. To this 
end, the previous data are extensively re-analyzed. Fi- 
nally we will discuss the origin of the short range spin 
order related to the incommensurate spin structure. This 
structure is close to the commensurate structure with the 
periodicity of (2^3 x 2\/3). 

Recently the magnetic properties of 2D antiferromag- 
net on the triangular lattice, NiGa2S4 {S — I for Ni^+) 
have been extensively studied .'^^^^^'^ There are several 
features which is common to stage-2 MnCl2 GIC, because 
of the same universality class, (i) The incommensurate 
short range spin order with nanoscale correlation is ob- 
served in the temperature range between 1.5 K and 15 
K. (ii) The spin glass phase appears below Tsg 10 
K). (iii) The magnetic specific heat shows no anomaly 
at Tsg- It shows double broad peaks at 13 and 80 K. 
(iv) The magnetic specific heat is proportional to in 
the low temperature side, (v) The entropy due to the 
broad peak around 13 K is 1/3 of the total entropy. It is 
surprising that in spite of difference in the detail of mag- 
netic and structural parameters, the magnetic behavior 
of stage-2 MnCl2 GIC is very similar to that of NiGa2S4. 
This may imply that these common features are a key to 
understanding the physics of 2D XY-like or Heisenberg- 
like antiferromagnet on the triangular lattice. 



stage-2 MnCI^ GIC 




FIG. 1: (Color online)Plot of specific heat C in the units of 
J/(Mn mol K) vs T for stage-2 MnCIz GIC. H ^ 0. Note 
that two kinds of data which were separately taken under the 
same condition, are shown. 



II. EXPERIMENTAL PROCEDURE 

The sample of stage-2 MnCl2 GIC was prepared from 
kish graphite by vapor reaction of the powdered MnCl2 in 
a chlorine-gas atmosphere of 740 Torr. The reaction was 
continued at 500°C for one month. The stoichiometry of 
the sample was determined from the weight-uptake mea- 
surement as Ci2.7oMnCl2. The specific heat was mea- 
sured by the thermal-relaxation method at low temper- 
atures down to 0.4 K in magnetic fields parallel to the c 
plane up to 20 kOe on field cooling, using a commercial 
apparatus (PPMS, Quantum Design). 



III. RESULT 

A. Determination of nonmagnetic specific heat at 

H = 

Figure [l] shows the T (temperature) dependence of spe- 
cific heat (C) in zero field for a stage-2 MnCl2 GIC based 
on single-crystal kish graphite, where the specific heat is 
expressed in the unit of J/ (Mn mol K). The specific heat 
increases with increasing T at low temperatures below 
5 K, becoming to flatten out between 5 K and 20 K, 
and it increases with further increasing T. There is no 
anomaly in C vs T around the spin freezing temperature 
Tsg (= 1-1 K), where the ZFC magnetization exhibits a 
relatively sharp peak. 

In Figs. [2j^a) and (b), we make a plot of C/T vs at 
H = for stage-2 MnCl2 GIC. It is found from Fig. [2] 



that C/T drastically decreases with increasing at low 
temperatures below 20 K. It shows a minimum around 
20 K and in turn it increases with further increasing T^. 

It is expected that the specific heat of stage-2 MnCl2 
GIC at low temperatures, C, is described by. 



C — Cman + Ci 



noirtag j 



with 



C„ 



Ce + Cg^jT + pT^ 



(1) 



(2) 



where Ce (=7^) and Cg (= (3T^) represent the contribu- 
tions of the conduction electrons and of the lattice to the 
specific heat, respectively, and Cmag is the magnetic con- 
tribution. Here we assume that the Deby temperature of 
stagc-2 MnCl2 GIC with stoichiometry Ci2.7oMnCl2 is on 
the same order as that of the pristine graphite. Therefore 
the T^ term of the Cnonmag may be negligibly small for 
T < 50 K. Then C/T is written as 



C 
T 



Cm a q , Cn omaq 



C 



T 



T 



mag 



+ j + f3T^ 



(3) 



In the temperature range where the magnetic contribu- 
tion Crnag/T IS negligibly small, the curve of C/T vs T^ 
is a straight line with the y-intercept 7 and the slope /3. 

In Fig. ^a.), we note that the data of C/T vs for 
20 < T < 50 K do not fall on a straight line, but is up- 
ward convex for any possible straight line corresponding 
to the plot of Cnonmag /T VS T^ . This may imply that 
the magnetic specific heat might exhibit a small broad 
peak centered around 40 — 50 K. In fact the existence of 
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FIG. 2: (Color onlme)(a) C/T vs for stage-2 MnCla GIG. 
H ^ 0. (b) Overview of the plot of C/T vs for stage-2 
MnGb GIG. H = 0. The solid line indicates the C„onmag/T 
vs r^, which is obtained from the procedure described in the 
text. 



the small broad peak has been confirmed from the AC 
calorimetry measurement by Koga and Suzuki.^ Such a 
shght deviation of C/T vs T^ for 20 < T < 50 K from 
a straight line makes it difficult for one to determine ex- 



actly the T dependence of Cnonmg ■ Here we assume that 
the straight line for the plot Cnomag/T vs T^ passes two 
points at r = 22.8 K and 50.5 K. It is expressed by 



a 



nonmag 

~T 



(4) 



with the parameters 7 and (5 given by 7 = 0.38807 (J/Mn 
mol K^), and (3 = 8.4397 x lO'^ (J/Mn mol K"*). The 
advantage of the choice of the straight line for Cnonmg /T 
vs T^ is that the plot of C/T vs is always above 



the straight line for Cnonmg /T vs T^. In other words, 
Cmag IS always positive. The disadvantage is that the T 



dependence of Cmag 
to the choice of the straight line 



well above 20 K is rather sensitive 
Fortunately, the T 



dependence of Cmag well below 20 K, is not so sensitive 
to the choice of straight line, since the magnitude of Cmag 
is much larger than that of Cnonmag- 



B. Magnetic specific heat Cmag vs T 



In Fig. ^^&) we show the T dependence of Cmag for 
-ff = at low temnperatures, from the subtraction 
Cnonmag from C. Since Cnonmag IS vcry Small at low tem- 
peratures, the uncertainty of Cmag is negligibly small. 
The magnetic specific heat Cmag shows a broad peak 
around 5 K. No anomaly in Cmag vs T is observed at 
the spin freezing temperature Tqq (= 1.1 K). Above T^p 
(~ 5 K), Cmag decreases as 1/T with increasing T. In 
Fig. [3]^b) we show the T dependence of Cmag between 
0.4 and 50 K. The magnetic specific heat exhibits a lo- 
cal minimum around 22 K, and a local maximum around 
41 K. The peak height of Cmag around 41 K is depen- 
dent on the choice of Cnonmag as a background, since the 
contribution of Cnonmag to the specific heat C is large. 
The anomaly at 41 K is closely related to the change of 
the spin symmetry from isotropic to XY like. Wiesler 
et al.t^ have reported the T dependence of the DC mag- 
netic susceptibility (x_l along the c plane and x\\ along 
the c axis). Both x_l and x\\ decrease monotonically 
with increasing T and exhibits no anomaly around 50 K. 
Although x_L > X|| at all temperatures, the anisotropy 
= X-L — X\\ becomes appreciable only below 50 K. 
This may indicate the onset of short-range spin correla- 
tion in MnCl2 layer. 



T dependence of magnetic specific heat Cmag 
the presence of H 



Figure |4] shows the T dependence of the magnetic spe- 
cific heat Cmag, where an external magnetic field H is 
applied along the c plane. The external magnetic field 
H is changed as a parameter; H = 0, 2, 10, and 20 kOe. 



It is found that C, 



mag IS strongly dependent on the mag- 
when the 



IV C 



netic field. As will be discussed in Sec 
magnetic field is applied to our system at T below Tsg, 
the spin glass phase is destroyed and changes into the 
paramagnetic phase above the critical field HdT). This 
is called the de Almeida-Thouless (AT) transition .I^^ The 
critical field Hc{T = 0) (defined as Hat) is evaluated as 
5.87 ± 0.47 kOe. This implies that the system is at least 
in the spin glass phase at low fields (iJ = 0, and 2 kOe). 
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FIG. 3: (Color online)Plot of Cnag vs T for T < 15 K. Cmag 
shows a broad peak well above Tsg (= 1-1 K). Cmag exhibits 
no anomaly at T = Tsg, but exhibits a broad peak around 
Tcp (— 5 K). (b) Overview of the curve of Cmag vs T for 
< T < 50 K. Cmag shows double broad peaks around 5 and 
41 K. The height of the peak around 41 K may be dependent 



on the choice of C„o 



as a background. Since C„o 



very small at low temperatures below Tsa, the T dependence 
of Cmag is independent of the choice of Cnonmag- 



D. T dependence of Cmag below Tsg 
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FIG. 4: (Color online) Magnetic specific heat Cmag vs T at 
low temperatures for stage-2 MnCb GIC. An external mag- 
netic field H is applied along the c plane (perpendicular to 
the c axis). H = 0,2, 10, and 20 kOe. 
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FIG. 5; (Color online) Plot of Cmag/T as a function of T 
where H is changed as a parameter. H = Q, 2, 10, and 20 
kOe. H is applied along the c plane. 



Figure Is] shows the plot of Cmag/T vs T for stage-2 
MnCl2 GIC at low temperatures below 4 K, where H is 
changed as a parameter. For H = and 2 kOe, the curve 
of Cmag/T VS T drastically decreases with increasing T. 
For H = 10 kOe, it gradually decreases with increasing 
T. For H — 20 kOe, in contrast, the curve of Cmag/T vs 



T starts to increase linearly with increasing T and shows 
a broad peak around T = 2 K. 

Here we examine the T dependence of Cmag below Tsg 
for each H. For H — 20 kOe, as shown in Fig. |6] the 
least-squares fit of the data of Cmag/T vs T for 0.45 K 
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FIG. 6: (Color online) Plot of Cmag/T vs T ioi H = 20 kOe. 
Cmag/T is proportional to T below Tsg- The solid line is a 
least-squares fitting curve. 

<T<TsG to 

Cmag/T = ttQ + ^qT, (5) 

yields the parameters 

ao = 0.0186 ± 0.0088 J/Mn mol 

and 

(5o = 0.6253 ± 0.0116 J/Mn mol 

This implies that the magnetic specific heat at iJ = 20 
kOe is a sum of linear T dependence and the depen- 
dence below Tsg- 

For i7 = 10 kOe, as shown in Fig.|7] Cmag is almost lin- 
early proportional to T at low temperatures. The least- 
squares fit of the data of Cmag vs T at = 10 kOe in 
the temperature range between 0.45 and 1.5 K to 

Cmag = «! + 5iT, (6) 

yields the parameters 

ai = 0.077337 ± 0.001986 J/Mn mol K 

and 

5i = 0.80474 ±0.0019118 J/Mn mol 

We note that such a Hnear T dependence of Cmag at 
low temperatures is also observed in typical spin glasses 
such as Cuo.QssMno.oiJ^ and Euo.4Sro.6S.'^ There are 
several models to explain the linear T dependence in spin 
glass. These models are discussed in detail by Fisher and 



Stage-2 MnClj GIC 

1 .4 I — \ — \ — I — I — I — I — \ — \ — I — I — I — I — I — \ — \ — I — I — 1— 




4 I I 

0.4 0.6 0.8 1 1.2 1.4 1.6 



T(K) 

FIG. 7: (Color online) Plot of Cmag vs T for = 10 kOe. 
The solid line is a bets-fitting curve (straight line). 

Herz.l^ A dependence of Cmag at low temperatures 
may indicates the presence of gapless and linearly dis- 
persive modes in two dimensions. Such dependnce of 
Cmag at low tempcraturcs is observed in NiGa2S4.'Sl 

For H = 0, as shown in Fig. [5| the upturn of the curve 
for Cmag/T VS T ioi H = Oe is clearly seen in the low 
temperature ranges. Surprisingly, Wang and SwensoiJ^ 
have predicted a peculiar property of the asymptotic low 
temperature specific heat, 

Cmag ^ exp(-^), (7) 

or 

HCmagT^)=HA)~^, (8) 

using a replica Monte Carlo simulation method and 
finite-size transfer-matrix method for the 2D ±J Ising 
spin glass model, where A is constant and A is an en- 
ergy gap between the ground state and the first excited 
state in spin glass phase. The validity of this form has 
been confirmed by Saul and Kardar,!^ Lukic et al.jl^and 
Katzbraber et al. ^ in spite of the difference in the phys- 
ical interpretation of A. The plot of \ii{CmagT^) vs 1/T 
should be a straight line with a slope A. 

Figures [s] shows the plot oi\n{CmagT'^) vs \/T below 
Tsg for stage-2 MnCli GIC ai H = Q. It is found that 
the curve of \Yi{CmagT^) vs 1/T becomes straight line for 
0.45<T<0.7K, as shown in Fig|8] The least-squares 
fit of these data to Eq. (|8| yields the parameters 

A = 1.41 ± 0.05 K, \tl{A) = 1.09 ± 0.06. 

The energy gap may disappear when the Zeeman energy 
{gj.lJ'BSH) as result of the apphcation of H to that sys- 
tem, becomes larger than the energy gap fc^A, where S 
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FIG. 8: (Color online) Plot of \n{CmagT^) vs l/T at low 
temperatures for stage-2 MnCl2 GIC. The straight line below 
TsG is a least-squares fitting curve denoted by Eq.Q with 
In A = 1.0928 and A = 1.407 K. 
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FIG. 9: (Color online) Entropy (= ^(r) - S{Tmin) vs T, 
which is calculated from the data of Cmax vs T. H = 0, 2 
kOe, 10 kOe, and 10 kOe. T^in = 0.45 K. 



= 5/2 and g± — 1.977. The critical magnetic field Hq is 
evaluated as 



Ho 



4.25 kOe. 



The value of Hq is on the same order as the critical 
magnetic field at T = K for the AT transitionp' 
iJ^T = 5.87 ±0.47 kOe. 



E. Entropy vs T 

Figure |9] shows the plot of entropy AS = S{T) - 
S{Tmin) as a function of T for stage-2 MnCl2 GIC for 

= 0, 2, 10, 20 kOe, where Tmin = 0.45 K. The entropy 
AS* is defined as 



A5 = 5(T)-5(r.„i„)= / 



dT'. 



(9) 



Since S — 5/2, the total entropy Stotai is equal to 



5, 



total 



i?ln(2S'+ 1) = i?ln6 = 14.898 J/Mn mol K, 



where R is the gas constant and R — 8.314472 J/Mn mol 
K. As shown in Fig. [9] the curve of AS* vs T at = 
is upward convex, while the curve of AS* vs T at H 
= 20 k Oe is downward convex. The curve of AS* vs T 
at = is almost straight line, reflecting the linear T 
dependence of Cmag at 7? = 10 kOe. Figure 10 shows the 
T dependence of AS* at = in the wide temperature 
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FIG. 10: (Color online) Overview of entyropy (= ^(r) — 
S{Tmin)) vs T, which is calculated from the data of Cmax vs 
T. H = 0. Tmin = 0.45 K. The entropy due to the broad 
peak around 5 K is 1/3 of the total entropy. 



region up to 50 K. The entropy AS" becomes flat between 
20 and 30 K. The value of AS" at 20 K is nearly equal 
to 1/3 of i?ln6. This factor 1/3 is also observed in the 
entropy in NiGa2S4,^* which may be feature common 
to the 2D AFT. There are two step-like changes around 
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may be a residual entropy at T = K because of highly 
frustrated nature of the system. This residual entropy is 
calculated as 



^residual 



Jtotal 



S{T 50 K) = 9.43 J/Mn mol K 



which corresponds to 63 % of the total entropy. The 
residual entropy occurs when the system can exist in 
a large number of ground states with the same zero- 
point energjf. Note that for highly frustrated pyrochlore, 



Dy2Ti207j251 the residual entropy is 1/3 of the total en- 
tropy. Here we note that even the residual entropy is 
finite, still the magnetic specific heat is zero at T = 
Kpfil since 



Sj-e 



idual 



lim S ■ 



TS 
lim 

T-*o T 



lim (T — 

T^o- dT 



S) 



lim C 



dS 

1 m T— = 0, 
T^o dT 



(12) 



FIG. 11: (Color online) Normalized entropy given by S/Smax 
vs X (= T/A), where S/Smax = f{x) defined by Eq.([ll]). A 
= 1.41 K. 



5 and 41 K corresponding to the peak temperatures of 
Cjnag VS T. It is Considered that there is small amount 
of entropy contributing to the total entropy above 50 K, 
because of the absence of anomaly in Cmag above 50 K. 

We now evaluate the entropy S(T) at low tempera- 
tures. To this end, we assume that Cmag is given by 
Eq.([7]). Then the entropy is calculated as 



S = 



c 



dT = 



:^exp(- — )rfr 



A 
TA 



2(r + A)exp(--), 



A. 
T' 



(10) 



leading to the scaling form of the entropy given by 

S fi T (l + a;)exp(-i) 

-fi-=^) = (11) 



max 



with X = T/A and S^ax = A/A'^ [=1.507 J/(Mn mol 
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K)], where f{x) is a scaling function of x. In Fig 
we show the plot of the normalized entropy S/Smax as 
a function of x {— T/A). As is clearly shown, the curve 
of S/Smax vs X IS Upward convex, which is similar to the 
the experimental curve of vs T at = (see Fig. [9]). 

When Tmzn = 0.45 K and A = 1.41 K, S{Tmin) can 
be evaluated as S{Tmin) = 1.507/(x = r™„/A) = 0.272 
J/Mn mol K from Fig. 11 Then the observed entropy at 
50 K is calculated as 

S{T = 50 K) = S{Tm,n) +AS = 5.47 J/Mn mol K 

This value of S{T = 50K) is much smaller than the total 
entropy (= 14.898 J/Mn mol K). The missing entropy 



IV. SUMMARY OF RESULTS FROM 
PREVIOUS WORKS AND FURTHER ANALYSIS 

A. Magnetic neutron scattering 

The MnCl2 layer forms a triangular Mn^+ lattice in 
stage-2 MnCl2-GIC. Primitive lattice vectors a and b 
have length |a| = |b| = 3.692 ± 0.0051, the same as 
for pristine MnCl2 (3. 693 A). The in-plane reciprocal lat- 
tice vector is given by a* and b* where |a*| = |b*| = 
47r/(A/3a) = 1.965A^^. The magnetic neutron scatter- 
ing results of stage-2 Mn Cl2 GI G have been reported 
previously by Wiesler et al! ^*^* ^°^ The in-plane magnetic 
scattering at 0.43 K is obtained by subtracting the inten- 
sity at 14.85 K from the corresponding intensity at 0.43 
K. Magnetic peaks are observed at |Q_l| = 0.522, 1.536, 
2.05, and 2.44 A'^ (see Figs. 5 and 6 of Ref.^Hl), where 
Q_L is the in-plane wave vector. These wave vectors are 
consistent with the model shown in Fig. [T2j According 
to this model, the magnetic Bragg reflections should ap- 
pear at the in-plane wave vectors Qj^ = G(h,k) ± ki, 
G{h, k) ± k2, or G(/i, k) ± (ki - k2), where G(/i, k) = 
ha* + kh* is the in-plane reciprocal lattice vector of 
MnGl2 layer, h and k are integers, and k^ {i = 1, 2) 
are the reciprocal lattice vectors of magnetic superlattice. 
Assigning |ki| = 0.522A-1 and |a*-ki| = 1.536A-\ 
the angle 9 between ki and a* can be estimated as 



cos 9 



|ki| 



kil 



2|ki||a* 



(13) 



or = 30 ± 2° . The magnetic Bragg reflections are pre- 
dicted to occur at |Q_l| = |a*-Hk2| = 2.03l^\ and 
|a* + ki| = 2.43lA~^, in good agreement with the exper- 
imental values, |Q_l| — 2.05 and 2.44 A~^, respectively. 

The temperature dependence of the scattering has 
been also studied by a series of in-plane scans at temper- 
atures between 63 mK and 7.5 K. Throughout the entire 
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FIG. 12: (Color online) Magnetic Bragg reflections (denoted 
by solid circles) in the in-plane reciprocal lattice vector. 
|a*| = |b*| = 1.9651"^ |ki| = Ikal = 0.522l-i. The K 
point is denoted by (a* +b*)/3 or 120° spin structure or \/3 x 
spin structure; |Qcom| = V3\a*\ /3 = 1.1345 A" \ The 
points denoted by open circles correspond to the commensu- 
rate structure 2^3 X 2^3 spin structure. These points are 
expressed by (a* + b*)/6; |Qcom| = ^3 |a*| /6 = 0.56721"^ 
Note that the location of magnetic reflections is close to, but 
not equal to the magnetic Bragg reflections corresponding to 
the 2%/3 X 2\/3 commensurate spin structure. The reciprocal 
lattice of the graphite layer is not included in this flgure. 



temperature range, the intensities fit well to a Lorentzian 
peak described by 



/ = 



InK 



(14) 



around |Q_l| ~ |ki|, where k is the inverse in-plane 
spin correlation length, t is the peak position, and Iq 
is the integrated intensity. In Fig. 13 we show the T 
dependence of the peak position r. The peak position 
T is nearly constant at 0.522A~^ below 0.45 K and in- 
creases rapidly with increasing T in the vicinity of Tsg 
(= 1.1 K). Above Tsg, t plateaus briefly at 0.5321^1, 
rising again above 4.5 K (close to the peak tempera- 
ture of C„iag vs T) toward the commensurate po- 
sition at IQcomI = \/3|a*|/6 = 0.5672l-i (see the 
original data of t vs T in Fig. 7(a) in Ref.^°), where 
Qcom = (a*+b*)/6. 

In Fig. [14] we show the T dependence of the inverse 
in-plane spin correlation length /c (see also the original 
data of K vs T in Fig. 7(c) in Ref.'l^. The value of n is 
0.0566A~^ at 0.45 K, corresponding to an in-plane spin 
correlation length ^ (= 1/k) of only ISA. The value of k 
slowly increases with increasing T. It does not show any 




FIG. 13: (Color online) T dependence of the center of the 
Lorenzian peak, r, for stage-2 MnCl2 GIC. The original data 
of r vs T were taken by Wiesler et al.'^ r — \Qcom\ = 
\/3 |t[*| /6 — 0.5672^4"^ corresponds to the commensurate in- 
plane spin structure with the periodicity 2\/3 x 2\/3. The 
solid line is a guide to the eyes. Tsc ~ 1-1 K- Tcp ~ 5 K. 



0.14 



Stage-2 MnCI GIC 



< 




FIG. 14: (Color online) Inverse in-plane correlation length k 
(= 1/^) vs T. ^ is the in-plane spin correlation length. The 
instrumental resolution limit for k is O.OieA"^ The original 
data of K vs r were taken by Wiesler et al.'^ The solid line is 
a guide to the eyes. Tsg ~ 1.1 K. Tcp ~ 5 K. 
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stage-2 MnCI GIC 




FIG. 15: (Color online) T dependence of Mzfc an Mfc for 
stage-2 MnCl2 GIC. See the detail of ZFC and PC cooling 
protocol in the text. The magnetization is measured &t H = 
100 mOe. The original data were taken by Matsuura et al.l^ 




FIG. 16: (Color online) Plot of the peak temperature (x' vs 
T in the presence of H) as a function of H for stage-2 MnCl2 
GIC. jff II c plane. The original data were taken by Kimi shima 
et al.'^ The sohd line is a best-fitting curve to Eq.( 15 1. 



anomaly at Tsg ■ It undergoes a step-like change around 
4.5 K, at which Cmag shows a broad peak. The value of 
K is 0.124A~^ at 4.5 K, corresponding to an in-plane spin 
correlation length ^ of only 8A. 



summary, the stage-2 MnCl2 GIC undergoes a spin glass 
phase transition at T = Tsg ■ 



AC magnetic susceptibility, de 
Almeida-Thouless line 



B. ZFC and FC magnetization 

The T dependence of zero-field cooled (ZFC) magne- 
tization and field-cooled (FC) magnetization for stage-2 
MnCl2 GIC has been reported by Matsuura et al.^ Their 
result is reproduced in Fig.[l5] In the ZFC cooling proto- 
col, the system is quenched from high temperatures well 
above TgQ down to 0.5 K in the absence of an external 
magnetic field. Then the external magnetic field {H — 
100 mOe) is applied to the system at T = 0.5 K. The 
ZFC magnetization is measured with increasing T from 
0.5 K to temperatures well above TgQ under the presence 
of H . The FC magnetization is measured with decreas- 
ing T from high temperatures well above Tsg down to 
0.5 K, in the presence of H (= 100 mOe). The T de- 
pendence of Mzfc a.nd MpG is very similar to that of 
typical spin glases. Mzfc shows a peak around T = 0.9 
K. Mzfc starts to deviate from Mfc below Tsg (~ 1-0 
K). In other words, the irreversible effect of magnetiza- 
tion occurs below Tsg- The value of Tsc thus determined 
from Fig. [15] is a little different from that of the AC mag- 
netic susceptibility reported by Kimishima et aLl^Here 
we define that Tsg is equal to 1.1 K. The magnetiza- 
tion Mfc bends, becoming almost flat just below Tjg, 
which is one of the significant features of spin glasses. In 



Kimishima et al.l^ have reported the T dependence of 
AC magnetic susceptibihty (the dispersion x') for stage 2 
MnCl2 GIC in the presence of an external magnetic field 
H along the c plane. The dispersion x' shows a peak at 
T = Tsg (= 1.1 K). This peak shifts to the low temper- 



ature side with increasing H. In Fig. 16 we show the H 
dependence of the peak temperature for stage-2 MnCl2 
GIC which has been previously reported by Kimishima 
et al.l^ It is predicted for typical spin glasses that the 
peak temperature Tsg{H) forms a de- Almeida-Thouless 
(AT)^^ line in the H-T magnetic phase diagram, which 
is defined by 



H — Hat[^ 



TscjH) 
TscH^Oy 



(15) 



where p is an exponent and p = 3/2, and Hat is the 
magnetic field at T = K. The least-squares-fit of the 
data of Tsg{H) vs H to Eq.(15l yields the parameters 



TsciH = 0) = 1.184 ± 0.013 K, 

p = 1.79 ± 0.10, and Hat = 5.87 ± 0.47kOe. 

The exponent p thus obtained is relatively close to p = 
3/2 for the AT line. This indicates that the low tem- 
perature phase below Tsc is a spin glass phase. Here 
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we note that the characteristic field Hat is theoretically 
predicted as (Katori and Ito)^^^ 



stage-2 MnCI^ GIC 



Hat 



ksTsaiH = 0) 



(m + l)(m + 2) gp-sS 



(16) 



where m ~ \ for the Ising symmetry and m — i for the 
Heisenberg symmetry, g is the Lande g-factor, fc^ is the 
Boltzmann constant, and /is is the Bohr magneton. The 
value of Hat is evaluated as 2.103 kOe for to = 1 and 
as 3.839 kOe for to = 3, where g = g± = 1.97 for stage- 
2 MnCl2 GIC. These values of Hat is smaller than tha 
derived from the analysis of the least-squares fitting. 



D. DC magnetic susceptibility 

The DC magnetic susceptibility of stage-2 MnCl2 GIC 
has been reported by Wiesler et alW The parameters of 
the magnetic properties are listed in Table|T] Both x\\ and 
X± (measured in the FC state) decrease monotonically 
with increasing T and exhibits no anomaly around 41 
K, where the magnetic specific heat Cmag shows a broad 
peak. X-L > X\\ at all temperatures, which suggests that 
spins align in the MnCl2 layer. The anisotropy = 
X± — X\\ becomes appreciable only below 50 K, a result 
with a shift of the g- values below 50 K, may indicate that 
onset of short-range spin order with Q = mode in the 
MnCl2 layer. 

As listed in Table |l] the Curie- Weiss temperature <d± 
along the c plane is equal to —9.1 K, and the Curie- Weiss 
temperature along the c axis is equal to —5.9 K. The 
negative sign of Q± and indicates that the intrapla- 
nar exchange interaction is antiferromagnetic. The peak 
temperature {Tcp ~ 5 K) of Cmag vs T is on the same 
order as |6||| and |6_l|. 



E. Line width of electron spin resonance 

The electron spin resonance (ESR) of stage-2 MnCl2 
GIC has been measured by Koga and Suzuki,^ and 
Suzuki et al.l^Fi gure|l7|shows the T dependence of ESR 
linewidth AH for stage-2 MnCl2 GIC when the external 
magnetic field H is applied along the c axis (|j). The line 
width AH is defined as the peak-to-peak width of the 
derivative ESR line shape. 

The T dependence of ATfy for if || c is consistent with 
that of Tx\\ for 6 < T < 20 K. This suggests that the 
Q = mode is dominant in these temperature ranges. 
As T decreases from 6 K (« |0||| = 5.9 K), the mode 
Q — Qincomm of the Spin fluctuations takes over the Q 
= mode, where Qmcomm is the wave vector correspond- 
ing to the in-plane incommensurate spin structure (close 
to the 2a/3 X 2^3 spin structure). The minimum of A77j[ 
around 6 K results from the competition between the Q 
= mode and Q = Qincomm mode. The magnitude 
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FIG. 17: (Color online) T dependence of the ESR line-width 
AH for stage-2 MnCl2 GIC. H \\ c. The ESR measurement 
was made by the conventional field-modulated X-band spec- 
trometer at frequency v = 9.42 GHz in the temperature range 
between 1.38 and 293 K. The original data were taken by Koga 
and Suzuki,'^ Suzuki et al.'^ 



(= \Qincomm\) decreases with decreasing T as is ob- 
served from the magnetic neutron measurements. The 
linewidth AH\\ diverges as T approaches Tsg from the 
high temperature side. 



V. DISCUSSION 

A. Nature of short range spin order in the vicinity 

of Tsg 

The origin of the incommensurate in-plane short range 
spin order is discussed, based on the model proposed 
by SakakibarEp^ who explains the magnetic structure 
of CeEu (donor-type stage-1 GIC), where CgEu has a 
hexagonal closed packed (hep) type structure. Here we 
assume that stage-2 MnCl2 GIC also have the hep type 
structure, which is schematically shown in Fig. |18[ For 
simplicity, only the positions of Mn^+ ions are shown. 
From the symmetry of the stacking sequence in hep- type 
structure, we assume that there are two kinds of stacking 
sequences such as ABABAB (Fig.[l8| and AC AC AC ■ ■ ■ 
in the weight of 1 to 1, where A, B, and C layers forms 
a triangular lattice. The positions of Mn^+ ions in the 
B and C layers are shifted by the translation vectors A 
and 2A with respect to those of atoms in the A layer, 
where 



A=^(2a + b). 
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TABLE I: Parameters of the structure and magnetism in stage-2 MnCk GICpEElIMlIIol 



parameter 


value 


description 


|a*| = |b*| 


1.965 


in-plane reciprocal lattice of MnCl2 layer 




0.5672 A"^ 


wave number for the 2\/3 x 2-\/3 spin structure 


\QlZl!n\ 


0.522 A"^ 


wave number for the incommensurate short range spin order at 0.45 K 


a 


3.692 A 


in-plane lattice constant of MnCb layer 


d 


12.80 A 


c-axis repeat distance 


Ox 


- 9.1 K 


Curie- Weiss temperature along the c plane 


©II 

II 


- 5.9 K 


Curie- Weiss temperature along the c axis 


Px 


5.831 


effective magnetic moment along the c plane 


P|| 


5.838 /^s 


effective magnetic moment along the c axis 


5± 


1.977 


g-factor along the c plane at room temperature 


511 


1.912 


g-factor along the c axis at room temperature 


Ji 


- 0.20 K 


antiferromagnetic n.n. intraplanar exchange interaction 


D 


0.97 K 


positive single ion anisotropy, showing XY symmetry 


TsG 


1.1 K 


spin freezing temperature 


S 


5/2 


spin of Mn^+ 




Mn(A) 



Mn (B) 



Mn(A) 



Mn (B) 



Mn(A) 



FIG. 18: (Color online) Simplified model for the hexagonal- 
closed packed (hep) type structure of stage-2 MnCb GIC, 
where only the positions of Mn^"*" ions are shown, a — 3.692A. 
d — I2.8OA. a, b and c are lattice vectors. 



The lattice vectors are described in terms of the Cartesian 
coordinates, as 




,c = (0,0, 2d) 



where a is the lattice constant {— 3.692A, d is the c-axis 
repeat distance {d — 12. 80 A), and c — 2d is the lattice 
constant along the c axis. The reciprocal lattice vectors 
are also described in terms of the Cartesian coordinates, 




;>0 ,b* = 



47r / 1 



V3a\ 2 ' 2 



,0 



2n 
2d 



(0,0,1). 



where a* = |a*| = An /{Via) and c* = |c*| = n/d . The 
K point in the Brillouin zone (see Fig. 12 1 is expressed 

by 



l(a* + b*)=g (1,0,0) 



a* -^(1,0,0). 



The reciprocal vectors a* and b* are denoted in Fig. 12 
where the K point is on the Qx axis. According to 
Sakakibara,^® the minimum energy of the system is given 

by 



= -NS^JiQ), 



with 



(17) 



(18) 



where Q is the characteristic wave vector, the N is total 
number of Mn^+ ions. Si is the classical spin vector of 
Mn^+ ion at the site, R^- = R^ — Rj, J(Ry) is the 
exchange interaction between Si at the Ri site and Sj 
at Rj site: J(Rij) = J(— Rjj). In our system, J(Q) is 
defined by 



J(Q) = Ji(Q) + MQ) + MQ) + J'(Q), 



(19) 



12 



with 

MQ) 
MQ) 



,aQx _„„/\/3 



2Ji[cos(aQ^) + 2cos(^^) cos( — aQy)], (20) 



3 

= 2J2[cos(\/3a(3y) + 2cos{-aQx) cos( — aQ^)], 

(21) 

= 2J3[cos(2a(52:) + 2cos(a(5a;) cos(-\/3a(5y)], (22) 
= 2J'cos(dQ,)[2cos("|^)cos(0|)+cos(^)], 

(23) 



where Ji is the nearest neighbor (n.n.) intraplanar ex- 
change interaction (antiferromagnetic in the present sys- 
tem), J2 is the next-nearest neighbor (n.n.n.) intrapnar 
exchange interaction, J3 is the next next-nearest neigh- 
bor (n.n. n.n.) intrapnar exchange interaction, and J' is 
the interplanar exchange interaction. The wavenumbers 
Qx, Qy, and Qz are defined as 



47r 



Ait 

h — a*h, Qy — k — a*k 



The in-plane spin structure is reaHzed as a result of the 
minimum energy state with a characteristic in-plane wave 
vector Q, as 



Si — S'cos(Q • Ri)e2; -I- 5'sin(Q • Hijey 



(24) 



on the triangular lattice site Ri, 

We need to find the location of the magnetic Bragg 
points in the Q space for each pair of Ji, J2, J3, and 
J'. The magnetic Bragg reflection appears at the vec- 
tor Q, where J(Q) has a maximum. To this end, we 
use a ContourPlot (equi-energy contour plot) of J(Q) in 
the Mathematica, where the points with the same J(Q) 
are connected, forming an contour line. When Ji = — 1 
(for simplicity the magnitude is normalized as 1, anti- 
ferromagnetic), J2 = J3 = J' = 0, the magnetic Bragg 
points appear at the K point, (l/3a* + l/3b*). The in- 
plane spin structure with Q at the K point is a 120° spin 
structure (a/S x \/2i) as is expected. 

We determine the magnetic phase diagram in the ( J2, 
J 4) plane where Ji — -1 and J3 = 0. The magnetic phase 
diagram consists of the K point reflection ( ^/?> x ^/?>) , the 
M point reflection (2x2), and (-\/7 x %/?) reflections, de- 
pending on the values of J2 and J4. The detail of the 
magnetic phase diagram will be reported elsewhere. In 
the limited case of Ji = -1, -0.02 < J2 < 0.08, and 
J' = 2, the magnetic Bragg points appear at the point 
(a*/6 -I- b*/6), or = ha* and Qy = k a* with h — 
V3/6 = 0.28675 and k = 0. In Figs, and (b), we 



0. In Figs. 

show the energy-contour plot of J(Q), where Ji = — 1, J2 
= , and J' — 2. It is clear that J(Q) has a maximum at 
h — 0.28675 and fc = 0, which is the Bragg peak from the 
commensurate in-plane spin structure. Note that the en- 



ergy contour plot of Fig. 19 'b) is similar to the schematic 
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FIG. 19: (Color online) Energy contour line of J(Q) for Ji = 
— 1, J2 = 0, Ji — 0, and J' = 2 in the {h, k) plane. Qx = a*h. 
Qy = a'k. = 0. h = 0.28675 and fc = for the \Qc ' - 
0.5672A-1. 



Ccomm 



Because of J' > 0, J' is a ferromagnetic interplanar in- 
teraction. The magnitude of J' is larger than that of Ji. 
This result seems to be inconsistent with our assumption 
that I J' I is much smaller than |Ji|. However, if the effec- 
tive interplanar exchange interaction J'^ff is used instead 
of J' and the magnitude of J'^j^j becomes large with de- 
creasing T, the above inconsistency may be ruled out. 
Below 30 K where the static magnetic susceptibility be- 
comes anisotropic, the short-range spin order gradually 
grows with decreasing T. This growth leads to the effec- 
tive interplanar exchange interaction between adjacent 
MnCl2 layers. 



(25) 



becomes 
Then 



with = 18^ at 0.45 K and a = 3.692^1. Even if 

I J' I is very small, the effective interaction 

large when become large with decreasing T 

J'^ff gives rise to a spin correlation between Mn^+ ions 

of adjacent MnCl2 layers. 

In the present stage, we cannot flnd any appropriate 
values of (J2, J4) for the magnetic Bragg peak at the 
incommensurate spin structure {h = 0.2656, fc = 0) cor- 
responding to IQincomml — 0.522A^^ from our method. 
Note that this situation may change when J3 is taken 
into account. 



diagram of the reciprocal lattice phase shown in Fig. 12 



Using Eq.(24l, the in-plane spin structure with a char- 
acteristic Q can be obtained on the triangular lattice 
without knowing the values of Ji, J2, J3, and J'. Fig- 
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FIG. 20: (Color online) (a) Commensurate in-plane spin 
structure with jQcommi ~ 0.5672A~^. 2\/3 x 2\/3 spin 
structure, (b) Incommensurate in-plane spin structure with 

IQincomml ~ 0.522.A 

ure |20[ a) and (b) shows the commensurate in-plane spin 
structures with |Qcomm| — 0.5672A^^, and incommensu- 
rate in-plane spin structure with |Qiricomm| = 0.522A^^. 

B. Nature of reentrant spin glass phase 

Below TsGi the short range in-plane spin order char- 
acterized by the incommensurate in-plane wavevector 
T — \Qincomm\ cocxists wlth the Spin glass phase. The 
in-plane correlation length slightly increases from 17 A at 
TsG to 18l 5a) at 0.45 K with decreasing T. The 
wavevector r = |Qmcomm| decreases from 0.535A~^ at 
TsG to 0.5221-1 gj- Q 45 ^ ^^^^ decreasing T, where 
IQcomml = 0.5672A-1. This means that the degree of 
the incommensurability of in-plane spin order with the 
Mn lattice is slightly enhanced well below Tsg- The in- 
plane spin correlation length almost remains unchanged 
below Tsg partly because of the highly geometrically 
frustrated-nature of the 2D AFT. 

The spin glass phase observed in our system is rather 
different from that of the dilute Ising antiferromagnet 
Feo.55Mgo.45Cl2. Wong et al.l^have shown that the spin 
glass behavior and the long-range antiferromagnetic (AF) 
order coexist in the low temperature phase, forming the 
reentrant spin glass phase (RSG). The system undergoes 



magnetic phase transitions at the Neel temperature Tjv = 
7.5 K and the reentrant spin glass transition temperature 
Trsg = 3.0 K. The AF phase with long range spin order 
exists between Tn and T^sg and the RSG phase exits 
below TjiSG- The AF order is long range (> 10^ A) both 
between T/v and Tj^sGi and below Te.sg- In other words, 
it is unaffected by the reentrant spin glass transition. 
The microscopic picture of the coexistence may consist of 
spins in an infinite AF network and spins frozen like spin 
glass. In our system, there are two kinds of spins, spins in 
the short range spin order and the frustrated spins. The 
spin direction of the frustrated spins are frozen below 
Tsg because of the spin frustration effect arising from 
the competing interactions. The short range spin order is 
not destroyed since this spin order is already established 
well above Tsg- 



VI. CONCLUSION 

The magnetic properties of stage-2 MnCl2 GIG are dis- 
cussed in terms of the experimental results of specific heat 
as well as ZFG and FG magnetization, AG magnetic sus- 
ceptibility, magnetic neutron scattering, and ESR. The 
short range spin order with associated with the incom- 
mensurate wave vector |Qmcomm| (= 0.5221-^ at 0.45K) 
appears below 5 K |6j_|, where the magnetic 

specific heat Cmag has a broad peak. The wavevec- 
tor IQmcomml IS cloSC tO the wavCVector (IQcomml = 

%/3 |a*| /6 = 0.56721-1) 

commensurate in-plane 
spin structure with (2\/3 x 2\/3) periodicity. The in- 
plane spin correlation length slightly increases with de- 
creasing T and reaches only 18l at 0.45 K below Tsg- 
In the low temperature phase below Tsg, the spin glass 
phase coexists with the short-range spin order associated 
with the incommensurate wave vector |Qmcomm|- The 
residual entropy at K is estimated to be 63 % of the 
total entropy because of highly frustrated nature of the 
system. The entropy due to the broad peak around 5 
K is 1/3 of the total entropy. The temperature {— 41 
K) at which Cmag shows a broad peak is the onset tem- 
perature of short-range spin order. These features are 
common to the systems with the same universality class, 
2D XY antiferromagnet on the triangular lattice. 
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